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Sakharov has proposed a suggestive model in which gravity is not a separately existing fundamen- 
tal force, but rather an induced effect associated with zero-point fluctuations (ZPF’s) of the vacuum, 
in much the same manner as the van der Waals and Casimir forces. In the spirit of this proposal we 
develop a point-particle-ZPF interaction model that accords with and fulfills this hypothesis. In 
the mode] gravitational mass and its associated gravitational effects are shown to derive in a fully 
self-consistent way from electromagnetic-ZPF-induced particle motion (Zitterbewegung). Because 
of its electromagnetic-ZPF underpinning, gravitational theory in this form constitutes an “already 


unified” theory. 


I. INTRODUCTION 


Gravitational theory, whether in its scalar Newtonian 
form or its tensor general-relativistic form, is recognized 
to be essentially phenomenological in nature. As such, it 


invites attempts at derivation from a more fundamental. 


set of underlying assumptions, and six such attempts are 
outlined in the standard reference book Gravitation, by 
Misner, Thorne, and Wheeler (MTW).! 

Of the six approaches presented in MTW, perhaps the 
most far-reaching in its implications for an underlying 
model is one due to Sakharov; namely, that gravitation is 
not a fundamental interaction at all, but rather an in- 
duced effect brought about by changes in the quantum- 
fluctuation energy of the vacuum when matter is 
present.”> In this view the attractive gravitational force 
is more akin to the induced van der Waals and Casimir 
forces, than to the fundamental Coulomb force. Al- 
though speculative when first introduced by Sakharov in 
1967, this hypothesis has led to a rich and ongoing litera- 
ture on quantum-fluctuation-induced gravity that contin- 
ues to be of interest. In this approach the presence of 
matter in the vacuum is taken to constitute a kind of set 
of boundaries as in a generalized Casimir effect, and the 
question of how quantum fluctuations of the vacuum un- 
der these circumstances can lead to an action and metric 
that reproduce Einstein gravity has been addressed from 
several viewpoints.* These are treated in some detail in a 
comprehensive review by Adler on gravity as a 
symmetry-breaking effect in quantum field theory.° 

On the basis of heuristic and dimensional arguments 
along general relativistic lines, Sakharov argues that in a 
vacuum-fluctuation model for gravity the Newtonian 
gravitational constant G should be determined by an ex- 
pression of the form 
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G ~—ao 
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where œ, corresponds to an effective Planck cutoff fre- 
quency of the vacuum zero-point-fluctuation (ZPF) spec- 
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trum.® In this approach, the small (but finite) value of the 
gravitational constant is an inverse reflection of the high 
(but not infinite) value of the high-frequency cutoff of the 
ZPF. 

In this paper we explore the Sakharov viewpoint on the 
basis of a conceptually simple, classical model (but in- 
cluding ZPF) in which matter, in the form of charged 
point particles (partons), interacts with the ZPF of the 
vacuum electromagnetic field. As part of this develop- 
ment the model predicts (1) to be precisely of the form 


5 1/2 
T c 
G = — — 0. = 
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In order to constitute a self-consistent, viable basis for 


rc’ 


4G (2) 


` gravitation, however, a first-order ZPF model for gravity 


must provide not only a basis for calculation of the gravi- 
tational constant G (shown to reflect the ZPF cutoff), but 
must also account for the genesis of the gravitational 
mass, and the attractive inverse-square-law force. In the 
particular version of the Sakharov hypothesis pursued 
here, the mass is shown to correspond to the kinetic ener- 
gy of ZPF-induced internal particle (parton) motion (ZPF 
“jitter,” or Zitterbewegung), while the force is found to 
be of a long-range retarded van der Waals type, associat- 
ed with the broad-spectrum ZPF radiation fields generat- 
ed by that same Zitterbewegung motion. 

To arrive at the above results, basically we simply as- 
semble together in a straightforward fashion previously 
published results regarding ZPF models of van der Waals 
and related effects in flat space-time.’ When this is done, 
one finds the leading term in the interaction potential, 
previously unexamined, to be Newton’s law with no free 
parameters to be fixed. ‘In such a fashion the 
identification of this term as gravitational emerges natu- 
rally from the concatenation of the previously published 
results. Yet further evidence for the correctness of this 
interpretation is provided in Sec. VI, where details of the 
pattern that emerges are presented in the context of a 
self-consistent coherent picture of the underlying dynam- 
ics of the gravitational force that conforms to the facts as 
we know them. 
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Il, VACUUM ZPF ELECTROMAGNETIC FIELDS (REF. 8) 


In the classical approach used here, point- 
particle-ZPF interactions are treated on the basis that 
charged point-mass particles interact with a background 
of random classical electromagnetic zero-point radiation 
with energy spectrum (as in the quantum-mechanical 
case) 


2 ha? 


T [dom 


ma ydo, (3) 
where the first factor in parentheses corresponds to the 
density of normal modes, the second to an average energy 
tw per mode. This treatment of quantum field-particle 
interactions on the basis of a classical ZPF constitutes an 
analysis technique known in the literature as stochastic 
electrodynamics (SED).? SED is a well-defined frame- 
work that has a long history of success in yielding precise 
quantitative agreement with full QED treatments of such 
topics as Casimir!°~ and van der Waals forces,’ topics 
directly related to the one pursued here. In this approach 
# appears in the above expression simply in the role of 
scale factor, without need of quantum interpretation, and 
all other appearances of ñ in the development can be 
traced back to its appearance in this expression. 

The spectral energy density represented by (3) formally 


p(o)do= | 2 


quantum field theory, which showed that, apart from an 
additional density-of-states factor [1+(a/ac)’], the spec- 
tral distribution seen by an accelerating observer assumes 
a thermal (Planck) form if one makes the identification 
T =ha/2mrck (k is Boltzmann’s constant; T, absolute 
temperature). In commenting on the additional density- 
of-states factor, Boyer points out that the additional -con- 


- tribution beyond the thermal (Planck) form is related to 


the space-time properties of an accelerating reference 
frame. This gives us a clue that, via the equivalence prin- 
ciple, this additional term can be related to the gravita- 
tional interaction.'® 

Of special interest here, therefore, is not the thermal 
term of interest in the original treatment, but rather the 
leading terms 


ho? + hoa? 
amc? = Awe? 
These indicate that an accelerated observer would see the 
background ZPF spectrum augmented by a term propor- 
tional to the square of the acceleration. Application of 
the principle of equivalence then indicates that the addi- 
tional spectral contribution seen in a frame with accelera- 


tion a should also be seen in a nonaccelerated frame with 
local gravitational field g produced by a gravitational 


p'(@)=p)(@) + Ap'(a)= (5) 


diverges as w°. It is generally assumed, however, that the 
spectrum is effectively cut off at a frequency roughly cor- 
responding to the Planck frequency, 


op=(c/AG) 2. 


As we shall see in the following sections, this assumption 
is supported by the line of development presented here. 

As a first step toward developing the hypothesized un- 
derlying ZPF basis of the gravitational interaction, we 
compare the forms of the spectral distribution of the ZPF 
of the electromagnetic fields as seen from unaccelerated 
and accelerated frames of reference. 

Of particular significance with regard to the spectral 
distribution in an unaccelerated frame, given by (3), is the 
fact of its Lorentz invariance, which derives specifically 


from the spectrum’s cubic dependence on frequency. The - 


cubic spectrum is unique in its property that delicate can- 
cellations of Doppler shifts with velocity boosts leaves 
the spectrum Lorentz invariant.! 

In an accelerated frame, on the other hand, the de- 
tailed balance of Doppler-shift cancellations is negated, 
with the result that the spectral distribution takes the 


form!’ 
o a 
pla)da= a Hlé | 
fio hia 
A 2  expl2rco/a)—1 ahs i 


where a is the proper acceleration relative to a Lorentz 
frame. This expression was originally derived by Boyer! 

on the basis of the random classical ZPF background as- 
sumed here. The purpose was to derive a result, first ob- 
tained by Davies'® and Unruh!” within the context of 


mass mg. Setting g=—a= —î, Gm, /r?, we obtain 
Fa Gm, 
Ap'(@)= (6) 
k Anes | r? 


Thus the principle of equivalence predicts an addition- 
al contribution to ZPF energy by gravitational mass, a re- 
quirement that must be met in any ZPF-based theory of 
gravitation. Since this additional contribytion of energy 
is electromagnetic in nature, we must ascribe to mass an 
appropriate electromagnetic-field-generating function, a 
point to which we return in Sec. IV. 


Ill. ZITTERBEWEGUNG MODEL 


--~As our basic point-particle~ZPF interaction model, we 


represent matter as a collection of charged point-mass 
particles (partons), in accordance with standard theory. 
In the development that follows it is not necessary to in- 
voke the details of particular parton representations (e.g., 
families of fractionally charged quarks) beyond certain 
general concepts, such as the “asymptotic freedom” of 
partons to respond to the high-frequency components of 
the ZPF spectrum as essentially free particles. It is 
necessary to focus to the charged-parton level, however, 
in order to represent properly the essentially equal mag- 
nitudes of proton and neutron contributions to gravita- 
tional mass. This accounts for the fact that charged and 
neutral matter participate equally in the gravitational in- 
teraction, based on underlying charged-parton interac- 
tions. 

We begin our discussion of particle-field interactions 
by examining the properties of a simple charged harmon- 
ic oscillator of natural frequency œg (corresponding to a 
binding force that is linear in displacement from equilib- 
rium), located at the origin and immersed in zero-point 


radiation. The (nonrelativistic) equation of motion for a 
particle of mass my and charge q, including radiation 
damping, is given by 


¥+qEzp . (7) 


mot+mosr= 3 
67E gC 


If we introduce the dipole moment, p=qr, and the damp- 
ing constant, 


T=q*/6reymoc? , 
we can write (7) in a form convenient for later discussion, 


p+osp=Tp't 6reoc T Ezp . (8) 
We will also be particularly interested in the kinetic ener- 
gy, 6=im oi, which becomes 

po, 

12rege T 
Once ‘written in this form, the oscillator equation of 
motion (8) and energy equation (9) refer only to the global 
properties of the oscillator (dipole moment p, natural fre- 
quency @p, and damping constant I), and do not involve 
individual mechanical properties of the oscillator such as 
charge or mass. 

With regard to the ZPF fields, the vacuum is assumed 
to be filled with a random classical zero-point elec- 
tromagnetic radiation whose Fourier composition under- 
lies the spectrum given in (3). Written as a sum over 
plane waves, the random radiation, which is homogene- 
ous, isotropic, and Lorentz invariant, can be expressed as 


1/2 
3 
Ezp( r,t) =Re Š f ake eae 


ieee +i0(k,o) ; (10) 
5 ee eee ee 


se eee, (9) 
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o=1 To 


x ek riot tidk,c) ; (11) 
where o=1,2 denote orthogonal polarizations, € and k 
are orthogonal unit vectors in the direction of the electric 
field polarization and wave propagation vectors, respec- 
tively, 0(k,o) are random phases distributed uniformly 
on the interval 0 to 27 ( independently distributed for 
each k,o), and œ= kce. 

It is at this point that we need to consider the 
correspondence between the above equations and the 
parton-ZPF interaction of interest. First, we treat the 
parton. as a two-dimensional (rather than three- 
dimensional) oscillator, drawing on previous studies that 
model spin as the “internal” angular momentum associat- 
ed with two-dimensional Zitterbewegung motion. 
Second, because we are interested primarily in the 
particle’s high-frequency Zitterbewegung response to the 
ZPF, whose spectral density increases as o, we may to 
first order neglect the binding-force term involving wo 
(asymptotic freedom as it relates to the ZPF). Finally, we 
also neglect the radiation-damping force in comparison 
to the inertial force and ZPF driving terms. (Alternative- 
ly, the restoring-force and radiation-damping terms can 
be carried through in the derivations that follow, and ap- 
propriate approximations introduced at the end, without 
change of result.) 

To obtain the expectation value of kinetic energy 
of the Zitterbewegung motion given by (9), (6) 
= (p?) /12regc°T, we follow a procedure due to Rue- 

0 Under the asymptotically-free-particle assumptions 
stated in the paragraph above, the x component of (8) 
takes the form (with % a unit vector in the x direction) 


dp, fo 1/2 
— = brege (Ezp), =6reoc I Re d’k(ER) a (12) 
dt o ZP 0 > J 3 Eme 
Integrating once with respect to time, we obtain 
2 1/2 
foes 3 T 32 (4.2 ikr iot +i6lk,o) 
È, = ree T Re J dt x Se k (2-2) Rares gikr—iot +i | 
2 1/2 ior __ , 3 A 
=6reeTRe| $ f dex) | [= gee) (13) 
o=1 8r €p mH 
The expectation value (p2) then follows from 
1/2 1/2 = 
fio fio’ eer] 
(p 2) =4Re(36nefeT? 3 d'k f PKT) Au 
4 2, z J J 8m €, BT Eo —io 
eiT 
x expli(k—k')-r +i6(k,o)—i6(k’,0")}) ; (14) 


where use of the complex conjugates and the notation + Re stems from the use of exponential notation. Equation (14) 


can, however, be simplified to 


(p2) =r T S f PeR _fo 1— coslor) 


eat Srey a? 


(15) 
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where averaging over random phases involves the use of 

( exp[i(k—k’)-r+i6(k,o)—iO(k',o’)]) =5,,,5(k—k’) . (16) 

With f d*k— f dO, J dk k?, and the angular integration in k taking the form 
2 
faar] > 2232 |= 
o=l 

we can rewrite (15) (with a change of variables to w=kc) as 


(p2) =12eořc r? i de w[1— coslwr)] 


= f da- RRP], (17) 


=6e fic Ta? [1+ 


(o,7P [1—cos(@,7)—(w,7)sin(a,7)] | , (18) 


where œ, is the assumed cutoff frequency, to be deter- originating in parton-motion response to the electromag- 
mined later.?! netic zero-point fluctuations of the vacuum. It is there- 

For @,T>>1, (p; 2Y reaches the value 7 ___fore simply a special case of the general proposition that 
aes 32a, i ooo - . -the internal kinetic energy of a system contributes to the 
(by) = beoc Tw (19) effective mass of that system.” This derivation of mass 


For the. two-dimensional. Zitterb ewegung motion as- 28 40 internal kinetic energy of motion is thus the first re- 


sumed, sult derived from the Zitterbewegung model. As will be 
ae i _... ..._._ Shown in later sections, it is this mass that is involved in 
(D) =2( p2) = etic Te? , (20) the gravitational interaction.” 
which, when substituted into (9), yields IV. ZITTERBEWEGUNG FIELDS 
(ey Thos l H =m (21) We turn our attention now to the fields generated by 


say, the x component of motion, we find that an assumed 
e ‘® time dependence substituted into (8) yields for the 
magnitude of any particular frequency component 


It is thus seen that the expectation value of the kinetic 
energy of parton Zitterbewegung motion reaches a finite 
magnitude, limited by the finite value of the (as yet un- 
determined) ZPF cutoff frequency. The energy calculated brege T x 
in this way is in the nature of the so-called “transverse B,(o)= ~ (8-2) Ezp(o) > (24) 
self-energy” (in QED) of a particle in response to the elec- i 
tromagnetic zero-point fluctuations of the vacuum.22 where the overtilde designates the magnitude of a fre- 
Since the energy associated with this Zitterbewegung quency component, and once again we have neglected the 
motion is an internal particle energy, that is, not directly binding and radiation-damping forces. This expression 
observable, we identify this energy as that corresponding can then be combined with the ZPF-field expressions (10) 


to the rest-mass energy of the particle, m, and (11), and the standard oscillating dipole formulas” 
46) Dior OO oaaae 
ME aS Daata e ena > eae E,(@)= Re = pe ~iIg (25) 
l 0 
As will be shown in Sec IV, @,=(1e5/4G )/ 2 in which 
case (22) reduces to “hæ : 
H,(@)= Re |—pe''F | , (26) 
ae? ea ie ease Sl Ss o 4r j 
G 


to yield expressions for the dipole fields generated by the 
In this view the particle mass m is of dynamical origin,  Zitterbewegung motion, viz., 


Sree’ r ho ie iot +ihk,o 
= STe T <2) 
E, z} P 5 f ak Brie, (€-X)e GJ, (27) 
2 67e,c°T 12 
— c 3 o ho 4.4), iot tiako) 
H, =— Re > dk (€-%)e 7? POR, ; (28) 
a dmr ai J o? BT Eo 


where, with f a unit vector in the direction j joining the dipole to the field evaluation point, 


1 
(kr? 


G=k3e* | (x2) xt |— 


+[32(f-X)—2] 


F=k RXR) |i 
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(29) 


(30) 


The energy density in the dipole-field distribution can be calculated from (27) and (28) as 


Awa = 460C Ez) + 4419( H3) 


2 
=pRe(> T 33 J ak f ae se 2 T žo |° | zw |? 
TEo | g=10=] o? 8m Eo Bm Eo 
X(E-R)(8": Rag, o')t +i0(k,o)—i0(k',o')] 
2 4 i 
+4Re( 5 a 5 5 fak f ak a T || 6reg T fo VII fo’ |” 
T | c=1o'=1 @ o”? 87o 8T Eo 
X(E-RNE-DIF| exp[—ilo—o')t +16(k,o)—i0(k',o a1) : (31) 
oe a a a a C Oa ae eee Pe 
This reduces to Since in the final analysis we are interested in the net 
3h T? contribution of a large collection of randomly oriented 


Awg= o p (FPH, (32) 


16r? o 


where we have made use of (17), and have averaged « over 
random phases by the use of 


( exp[ —il@— o )t +i6(k,o)—i6(k’,o’)]) 


=8 bou k-k’). (33) 


The part of the integrand in parentheses in (32) can be 
evaluated by means of (29) and (30), and yields 


\F/2+|G/? i 


=k 2 sin <a 4cos ne 4 cos ve sin 2y 


» (34) 
(kr)? (kr)* (kr) 


where Ņ% is the angle measured from the dipole-motion 
unit vector ¥ to the evaluation-point unit vector f. The 
first term proportional to 1/r? constitutes the radiation 
field associated with the ZPF-driven dipole. As shown 
previously by Boyer,”° this radiation just replaces that be- 
ing absorbed from the background, on a detailed-balance 
basis with regard to both frequency and angular distribu- 
tion, and therefore does not result in an incremental 
change to the background. Of the two remaining (induc- 
tion) field terms, the 1/r* term predominates at large dis- 
tances, and is therefore the one of interest here. Desig- 
nating the term of interest by a prime, we have 


A 3c T? cos*p f D, 
4q?r* 0 


This expression, obtained on the basis of considering a 
single (x) component of motion, must be doubled to take 
into account the contributions of the two (independent) 
degrees of freedom in the model. This leads then to an 
overall spectral density 


3fic lo cos" 


rr’ 


(35) 


ae 


individual particle motions, we average over the solid an- 
gle to obtain 


Acr? 
2mrt ` 


Since ee to (22) there is a relationship between IT 


Ap!) (o)= rear Apl2ar? siny d= (37) 


_ and the particle mass m for ZPF-driven Zitterbewegung 


motion, (37) can also be written 


co 


2ňotrt © 


Api (@)= (38) 


Zitterbewegung motion therefore leads to the generation 
of an electromagnetic field distribution in proximity to 
the mass that is proportional to frequency times mass 
squared, divided by r*, and upon detailed examination is 
found to be half electric, half magnetic. According to (6), 
moreover, a field of just this form is required by the prin- 
ciple of equivalence. Under the assumption that the 
gravitational and rest masses are identical (m,=m), (6) 
and (38) can be equated to obtain the cutoff frequency 


5 11⁄2 
ire? 
v= |e (39) 
which satisfies the Sakharov condition (1). In terms of 


the cutoff frequency ,, (39) can be inverted to yield the 
gravitational constant G in the form of the second (nonin- 
dependent) Sakharov condition, namely, 


(40) 


We see therefore that the principle of equivalence re- 
quires a certain modification of the ZPF background by 
gravitational mass, and that the Zitterbewegung model of 
mass implies a similar modification. Furthermore, the 
compatibility of the two conditions on a precise quantita- 


“tive basis requires only the equivalence of gravitational 
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and rest masses, and a limiting cutoff frequency for the 
ZPF background on the order of the Planck frequency. 
This specification of the cutoff frequency @,, and its rela- 
tionship to the gravitational constant G, is thus the 
second result derived from the Zitterbewegung model. 
Thus we have a derivation that yields the relationships 
postulated to exist 
breaking models of gravity as a symmetry-breaking 
effect.© 


V. GRAVITATIONAL FORCE, 


The derivations in Secs. III and IV dealt essentially 
with the characteristics of single masses. In this section 
we investigate the interaction between two such masses. 

The starting point is equations of the form (8), written 
for two masses, but modified to take into account the fact 
that each mass experiences not only the background ZPF 
field, but also the ZPF-driven dipole field of the other 
mass. The procedure followed here is precisely that 
developed by Boyer for the derivation of the retarded van 
der Waals forces at all distances between a pair of polar- 
izable particles.27. Therefore we need only outline the 
procedure as it applies to our case. : 

Two masses A and B (taken here to be equal for ease of 
discussion) are assumed to be located at positions r, and 


rg, respectively, with B located a distance R from A, 


along the positive z axis of a coordinate system centered 
at A. The equations of motion, generalized from (8), take 
the form. l 


Pi +arp 4 =p, +6rege T Ezplt 4t) 


+ 6reqc TE ypg(t 4t), (41) 
Ps tops =I pz + 6reg T Ezplrp,t) 
+6e€9¢*7TE,4(rpg,t) , (42) 


where E,,(r ,,¢) is the dipole electric field at the position 
of particle A, due to the motion of particle B, and so 
forth. 


As in the derivations of previous sections, an assumed 


e?! time dependence yields (for the magnitude of any 
particular frequency component) equations of the form __ 


in dynamical scale-invariance- 


H.E.PUTHOFF a 2 


we neglect the binding and radiation-damping terms 
[terms in œ% and F in (45)]. Furthermore, on the scale of 
interest in gravitation (distances large compared with the 
wavelengths of the predominant Zitterbewegung frequen- 
cies), in calculating the force we retain only the radiation 
field term œ 1/kR. With these assumptions, (45)—(47) be- 


come 


D™4=—o@ , , (48) 


3eikR 
i= = i Io? ER” ne4=0. 

Solutions for J 4,(@) and Bp,() are obtained straight- 
forwardly? from (43) and (44), subject to the conditions 
(48) and (49). It is then necessary to construct from these 
solutions the appropriate (van der Waals) forces between 
the pair of particles A and B. From classical theory the 
force on a dipole is given by 


(49) 


F=(p-VE+ 2 xB : (50) 


Following Boyer,” we note that the time-averaged 


force can be transformed as 


(F)=((p- IE) + (42 xB) 


=((p-VJE)+(px(VXE)), (51) 
where the transformation of the second term 
ap j= Tdp 
(4 xB)=a J, g xBdt (52) 


‘follows from an integration by parts and substitution 
from Maxwell’s equations, noting that the end-point 
terms do not contribute to the time average. Manipula- 
tion of the vector operators in (51) then leads to the rela- 
tively compact form 


(F)=(p, VE, +p, VE, +p,VE;) . (53) 


One further reduction is possible. The symmetry of 
the geometry dictates that the average force on a dipole is 
along the z axis joining the two particles. Therefore the 


expression for the average force finally simplifies to 


DP 4, (@)+7,5p,(@)=67€9¢*TE zp. (ar 4) » 


DB;,. ( (a) th Falo = bree TE gp, (a, Tz), (44) 


and similar equations for the y and z còmponénts, where l 


(43) 


AF)=(F,)= (0. a =) (54) 


a P az |? a 


- - -In determining the average force on, say, dipole B, we 


must keep in mind that the electric field at B consists of 
both the background ZPF field and the dipole field due to 
particle A, so that, for example, the first term in (54) be- 
comes 


ð 3 
ERAON) = (Ps. Š Erlan) 


(55) 


D =—0 toii ro , paaa E 
i 1 i 1 
ae fR |1 —_—- 46 
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The latter expressions for 7,, etc., are derived for the 
geometry under consideration from the dipole expres- 
sions (25) and (29). 

Consistent with the assumptions of previous sections, 


ð 
+ (Ps. Br edz (rp o) , 


and so forth. 
The mathematics of carrying out the averaging then 
proceeds term by term as in Ref. (27), using (as in previ- 


ous sections) the +Re and complex-conjugate notation, 


including averaging over random phases. The result for 
the special case of interest here (radiation field only, bind- 


ing and radiation-damping forces neglected) is, in terms. 


of the potential U,” 


3yp2 u —2uR . 
aat Re f" au? ; u=— 22, (56) 


U = 
4 r 0 R? 


The only difference here as compared to the derivation in 
Ref. 27 (aside from the specialization to the radiation 
field term) is the use of a finite cutoff frequency. 

This result is derived here and in Ref. 27 on the basis 
of point particles interacting with a classical zero-point 
field. For those who might be more familiar with stan- 
dard quantum calculations, this result has also been ob- 
tained by Renne” from quantum-electrodynamic calcula- 
tions using a nonrelativistic-quantum-oscillator model, 
and by Casimir and Polder’! using fourth-order perturba- 
tion theory in quantum electrodynamics. 

Equation (56) was derived for the case in which three 
degrees of freedom for particle motion are assumed. For 
the two-dimensional Zitterbewegung motion assumed in 
our case (N =2), geometrical considerations require that 
U be reduced by a factor (N/3)’=4 (see Appendix A). 
With this taken into account the solution to (56) becomes 
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y=- Line) te sath (57) 
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With the potential thus defined, the force is obtained 
from 


(59) 


We see therefore that the potential has the desired 1/R 
dependence required for gravity, modulated by a fine- 
structure overlay of the form [(sinR)/R] which has a 
spatial periodicity characteristic of the cutoff (Planck) 
frequency (~ 1078 cm). If we extract the leading (nonos- 
cillatory) term, we find for the potential and force 
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A careful examination of the details of averaging over 
the rapid spatial variation (see Appendix B) indicates that 
the particle experiences an average force (F} given by 
the leading term in (61). With I given by (23) and œ, by 
(39), (F } can then be written in Newton’s law form (with 
no adjustable parameters required), 


2 
A oe he EE A (62) 


39 GRAVITY AS A ZERO-POINT-FLUCTUATION FORCE 


2339 


This derivation of Newtons’s law, which expresses gravi- 
ty as a van der Waals force, is thus the third and final re- 
sult derived from the Zitterbewegung model. 


VI. DISCUSSION 


We begin our discussion by reiterating the logic flow of 
the approach pursued here. The basic thesis is the Sak- 
harov proposal that gravity is not a separately existing 
fundamental force, but rather a residuum force derived 
from zero-point fluctuations of other fields in the manner 
of the Casimir and van der Waals forces. Particularizing 
this hypothesis to the ZPF of the vacuum electromagnet- 
ic field, we identify the gravitational force as the van der 
Waals force associated with the long-range radiation 
fields (as opposed to the usual shorter-range induction 
fields) generated by the Zitterbewegung particle (parton) 
motion response to the ZPF of the electromagnetic field. 
The steps in this identification are three. 

First, particle mass is defined as the “internal” (that is, 
unobserved) kinetic energy of Zitterbewegung motion. Its 
value is set by a radiation damping constant I intrinsic to 
the particle, in conjunction with the value of the univer- 
sal cutoff frequency w.. Second, the value of the cutoff 
frequency œ, is determined by the equivalence principle. 
This principle sets a requirement that an expected addi- 
tional contribution to the free-space Lorentz-frame ZPF 
spectrum, viewed from an accelerated frame, is to be 
equated to a similar contribution expected in an unac- 
celerated frame, but in a gravitational field (that is, near a 
particle performing Zitterbewegung motion, and thereby 
generating a mass term). Third, a straightforward calcu- 
lation of the long-range van der Waals force associated 
with the radiation-field-correlated motions of such parti- 
cles, whose parameters are determined by the above two 
steps, leads to Newton’s law with no free parameters to 
be fixed. 

With a detailed theory in hand, certain attributes of 
the gravitational interaction become explicable in funda- 
mental terms. As mentioned earlier, the relative weak- 
ness of the gravitational force is due to the fact that the 
coupling constant determined by (40), G =2c>/ho?, 
reflects as the inverse square the high value of the ZPF 
cutoff frequency. With regard to the attractive nature of 
the force, this is simply a reflection of a property typical 
of van der Waals forces in general. The fact that gravita- 
tional interaction is characterized by a unipolar (single- 
valued) “charge” (mass) can be traced to a (positive only) 
kinetic energy basis for the mass parameter. 

The lack of shielding effects in gravity can also be 
comprehended on a rational basis. As understood here, 
this is a consequence of the fact that ZPF “noise” (quan- 
tum noise) in general cannot be shielded, a factor which 
in other contexts sets a lower limit on the detectability of 
electromagnetic signals. Specifically, in the case of ordi- 
nary electromagnetic shielding, macroscopic materials 
constitute dense boundaries that substantially alter field 
distributions, with shielding one consequence. In the 
gravitational case as modeled here, however, matter con- 
stitutes a dilute particle gas in an essentially high- 
frequency Hohlraum. As a result, particle-ZPF interac- 
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tions have negligible effect on the overall field distribu- 
tion (hence a lack of shielding), while nonetheless permit- 
ting particle-particle interactions that lead to an attrac- 
tive potential. 

In addition, implicit in the development pursued here 
are issues that extend beyond the gravitational interac- 
tion, such as mass renormalization and possible ZPF- 
induced contributions to the binding forces within the nu- 
cleus.2429 However, of the many inferences that can be 
drawn from this study, the most important is simply the 
fact that it is possible to carry through the basic Sak- 
harov program, namely, to uncover a basis for gravity in 
the ZPF of other (nongravitational) fields. In particular, 
we have been able to explicate a first-order model’? based 
on the ZPF of the vacuum electromagnetic field alone, 
once we take into account its effects on particle motion. 
The model thus details an electromagnetic basis for gravi- 
ty. Assuming the model is a proper representation of the 
gravitational interaction, the “already unified” aspect of 
the model would seem to mitigate against canonical at- 
tempts at unification of gravity as a separate force, or 
quantization of gravity as a separate field, in favor of a 
viewpoint more aligned with that presented here. 

It is therefore seen that a well-defined, precise quantita- 
tive argument can be made that gravity is a form of 
long-range van der Waals force associated with particle 
Zitterbewegung response to the zero-point fluctuations of 
the electromagnetic field. As such, the gravitational in- 
teraction takes its place alongside the short-range van der 
Waals forces and the Casimir force as related phenomena 
which emerge from the underlying dynamics of the in- 
teraction of particles with the zero-point fluctuations of 
the vacuum electromagnetic field. 
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APPENDIX A: DIMENSIONAL REDUCTION FACTOR 
FOR COUPLING CONSTANT 


As discussed in Sec. V, the attractive force between 
particles derives from coupling between ZPF-induced di- 
pole motions of the particles involved. The nature of the 
coupling is simplified somewhat in that the coupling 
takes place only between corresponding components of 
the motions of the two particles; that is, between J4, and 
Pex» Day and Pgp, and P 4z and Pg, This can be seen in 
the form taken by (43) and (44). 

In the derivation which leads to expression (56) for the 
interaction potential U, reference to the original deriva- 
tion in Ref. 27 shows that the potential is proportional to 
the quantity (n? + Ny 2 +72), where ny, etc., are as given 


here in (46) and (47). That derivation assumes full three- 


dimensional motion, with 4x coupling to Pp, By, to 
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For the two-dimensional Zitterbewegung motions po- 
sited here, however, we may for convenience analyze the 
couplings on the basis of assuming a random distribution 
of one-third each of x-oriented dipoles (y-z motion), y- 
oriented dipoles (x-z motion), and z-oriented dipoles (x-y 
motion). For two x-oriented dipoles, both y and z 
motions couple, yielding a contribution proportional to 
n +n. For an x-oriented dipole (y-z motion) coupling 
to a y-oriented dipole (x-z motion), only the z components 
couple, leading to a contribution proportional to 72; and 
so forth, For the nine possible dipole-pair combinations, 
straightforward enumeration of the possibilities leads to 
an average coupling factor proportional to 


(492 +492 +492) /9=4 2 +05 +077) « 


As a result, for the two-dimensional motions of interest 
here, a reduction factor of £ is to be applied to the value 
of the coupling constant obtained for the general three- 
dimensional case without constraints. 


APPENDIX B: AVERAGE FORCE 


In the Zitterbewegung model of gravity, the two- 
particle interaction potential based on the radiation field 
van der Waals effect is given by (57), repeated here, 

2 
sink 
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As seen, this expression can be factored into two parts, 
one with a slow spatial variation, 1/2, and one with a 
rapid spatial variation (on the order of the Planck wave- 


Jength), [(sinR)/R]?. Of interest in the gravitational in- 


teraction is not the rapidly varying component, but rath- 
er an average value, averaged over a distance large com- 
pared to the Planck wavelength. With the potential 
given by (B1), the (normalized) force is given by 
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As particle separation changes by an amount AR, the 
corresponding change in potential is given by 


R HAR 
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(B3) 
Assuming integration over a full cycle of the Planck vari- 
ation so that cos[2(R; +AR]= cos(2R;), and recogniz- 
ing that AR <<R; so that (R; ARP AR} + 3R7AR, we 
find that (B3) simplifies to 
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The change in potential, integrated over a cycle, is seen 
from (B4) to be sensitive to where in the cycle, 0; =2R;, 
the integration was begun. The average change in poten- 
tial is therefore determined by averaging over the range 
i | 


of possible initial starting points within the cycle, namely, 
an ER; Salan +1), where n is an integer n >>1. 

By reference to standard tables of integrals? we find, 
using (B4), 
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Substitution of the limits of integration, with the recognition that n >>1 implies that (n +1) =~n?+pn?”™!, then leads 
to 


3 5 
(ay) = YARAR onn T p en) |, (B6) 
(mn)? m 3! 5! 

But the term in large parentheses is recognized to be sin(2zn) =O, so that (AU) becomes 

(au) = LAE | (B7) 

R 
from which the average force can be calculated as 
ma WAU) De Stee 

(F) me Re oe (B8) 

The actual (unnormalized) force, F = —dU/0R, is recovered from the above with the aid of (58), yielding 
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